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Abstract 

We study the effect of nuclear equation of state on the r-mode instability of a rotating neutron 
star. We consider the case where the crust of the neutron star is perfectly rigid and we employ 
the related theory introduced by Lindblom et al. Q]. The gravitational and the viscous time 
scales, the critical angular velocity and the critical temperature are evaluated by employing a 
phenomenological nuclear model for the neutron star matter. The predicted equations of state for 
the /3-stable nuclear matter are parameterized by varying the slope L of the symmetry energy at 
saturation density on the interval 50 MeV < L < 110 MeV. The effects of the density dependence 
of the nuclear symmetry energy on r-mode instability properties are presented and analyzed. A 
comparison of theoretical predictions with observed neutron stars in low-mass x-ray binaries is also 
performed and analyzed. 
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I. INTRODUCTION 



Mi 



The oscillations and instabilities of relativistic stars (2|-|4j gained a lot of interest in the 
last decades because of the possible detection of their gravitational waves |5-14|. Especially 
neutron stars may suffer a number of instabilities which come in different flavors but they 
have a general feature in common, they can be directly associated with unstable modes of 
oscillation. In the present work we concentrate our study on the called r-mode instability. 
The r-modes are oscillations of rotating stars whose restoring force is the Coriolis force. The 
gravitational radiation-driven instability of these modes has been proposed as an explanation 



young neutron stars and of accreting 
1]. This instability can only occur when 



for the observed relatively low spin frequencies o 
neutron stars in low-mass X-ray binaries as well 
the gravitational-radiation driving time scale of the r-mode is shorter than the time scales 
of the various dissipation mechanisms that may occur in the interior of the neutron star. 

A very interesting problem is the consideration of the effect on r-mode instability due 
to the presence of a solid crust in an old neutron star. It is proved that the presence of a 
viscous boundary layer under the solid crust of a neutron star increases the viscous damping 
rate of the fluid r-modes Q, [l5]. Actually, the presence of a solid crust has a crucial effect 
on the r-mode motion and following the discussion of Andersson and Kokkotas this effect 
can be understood as follows: Based on the perfect fluid mode-calculations we anticipate the 
transverse motion associated with the mode at the crust-core boundary to be large. However, 
if the crust is assumed to be rigid the fluid motion must essentially fall off to zero at the 
base of the crust in order to satisfy a non-slip condition (in the rotating frame of reference). 



Firstly, Bildsten and Ushomirsky 



15] , found that the shear dissipation in the viscous 



boundary layer between the solid crust and the fluid core decreases dramatically the viscous 
damping time in cold old neutron stars as well as in hot young neutron stars. They concluded 
that the r-mode instability is unlikely to play a role in old, accreting neutron stars where 
in hot young neutron stars the boundary-layer damping mechanism limits the ability of the 
r-mode instability to reduce the angular momentum of the star and consequently to produce 



detectable amounts o: 
Anderson et al. 



gravitational radiation. 



16] used different neutron-star parameters to obtain significantly different 



results for the critical frequency of the onset of r-mode instability. They found the critical 
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velocity to be about 40% lower than the estimate of Bildsten and Ushomirsky and inferred 
that the r-mode instability is likely to be the mechanism that limits the low-mass x-ray 
binaries (LMXBs) spin periods and those of other millisecond pulsars as well flf] ]. 

Rieutord [17| improved the model of the boundary layer and found that the critical velocity 
agrees rather closely with the original estimates of Bildsten and Ushomirsky. Lindblom [lj 
improved previous estimates of the damping rate by including the effect of the Coriolis 
force on the boundary-layer eigenfunction, using more realistic neutron-star models. They 
concluded that if the crust is assumed to be perfectly rigid, the gravitational radiation 
driven instability in the r-modes is completely suppressed in neutron stars colder than about 
1.5 x 10 8 K and also found that the r-mode instability is responsible for limiting the spin 
periods of the LMXBs. 

Recently, Wen et al. 18] studied the sensitivity of the neutron star r-mode instability 



window to the density dependence of the nuclear symmetry energy. Employing a simple 
model of a neutron star with a perfectly rigid crust constructed with a set of crust and core 
equations of state that span the range of nuclear experimental uncertainly in the symmetry 
energy, they concluded that smaller values of the slope parameter L of the symmetry energy 



help stabilize neutron stars against runaway r-mode oscillations. Vidana [19] analyzed also 
the role of the symmetry energy slope parameter L on the r-mode instability by using both 
microscopic and phenomenological approaches of the nuclear equation of state. He showed 
that the r-mode instability region is smaller for those models which give larger values of L. 
Alford et al. [20J studied the viscous damping of r-modes of compact stars and analyzed the 
regions where small amplitude modes are unstable to the emission of gravitational radiation. 
They show that many aspects, such as the physically important minima of the instability 
boundary, are surprisingly insensitive to detailed microscopic properties of the considered 
form of matter. 

The motivation of the present work, in view of the above studies, is to study extensively 
the nuclear equation of state (EOS) effect on the r-mode instability of neutron stars with 



a perfectly rigid crust l|, ll8|, |2l|, |22j. Actually EOS affects the time scales associated with 



n 



the r-mode, in two different ways. Firstly, EOS defines the radial dependence of the mass 
density distribution p(r), which is the basic ingredient of the relevant integrals (see Sec. II). 
Secondly, it defines the core-crust transition density p c and also the core radius R c which is 
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the upper limit of the mentioned integrals. We employ a phenomenological model for the 
energy per baryon of the asymmetric nuclear matter having the advantage of an analytical 
form. By suitably choosing the parametrization of the model we obtain different forms from 
the density dependence of the nuclear symmetry energy by varying the slope parameter L 
on the interval 50 MeV < L < 110 MeV. The calculated EOS concerns the neutron star 
core (from the center of the star up to the crust-core interface). For the neutron star crust 



we employed the EOS taken from the previous work of Feynman, Metropolis and Teller 



23] 



and also from Baym, Pethick and Sutherland [24|. The transition pressure at the edge of 
the core is calculated by employing the thermodynamic method. It is found that in a good 
approximation the quantities connected with the r-mode instability are related directly with 
the slope parameter. The effects of the stiffness of EOS are studied on the variational and 
viscous time scales as well as on the critical angular velocity and critical temperature. Finally 
we investigate the case whether the observed LMXBs should constraint the slope parameter 
L. 

The article is organized as follows. In Sec II we review briefly the stability of the r-modes. 
Section III contains the employed nuclear physics model focusing on the presentation of the 
nuclear symmetry energy and the thermodynamic method applied for location of core-crust 
interface. The results are presented and discussed in Sec. IV and Sec. V summarizes the 
present study. 

II. STABILITY OF THE R-MODES 

The r-modes evolve with time dependence e lult - t / T as a consequence of ordinary hydrody- 
namics and the influence of the various dissipative processes. The real part of the frequency 
of these modes, u, is given by [3] 

(m — l)(m + 3) 



-a i 



m + 1 

where Q is the angular velocity of the unperturbed star. The imaginary part 1/r is de- 



termined by the effects of gravitational radiation, viscosity, etc. [l|, [25|, |26j. In the small 



amplitude limit, a mode is a driven, damped harmonic oscillator with an exponential damping 
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time scale 

1111 

V{Q/T) ~ 777> + T bv (V,T) + r v (Q,TY U 
where tqr, t^, t v are gravitational radiation, bulk viscosity and shear viscosity times scales. 
Gravitational radiation tends to drive the r-modes unstable, while viscosity suppresses the 



instability. More precisely dissipative effec 
(i.e., the mode is stable) as long as r > 



;s cause the mode to decay exponentially as e 



-t/r 



25]. 



The damping time t% for the individual mechanisms is defined in general by [l|, [ji], 

1 _ 1 (dE\ 
Ti = ~2E Vdt), 



- ■ (3) 



In Eq. ([3]) the total energy E of the r-mode is given by 

r-R 



1 P u / r \ 2m+2 

E = t ? ir ™-Wj o p(r)Q dr, 



(4) 



where a is the dimensionless amplitude of the mode, R and Q are the radius and the angular 
velocity of the neutron star respectively and p(r) is the radial dependence of the mass density 
of the neutron star. Similar expressions holds for the dissipation rate (dE/dt)i j^ . 

The damping time scale due to viscous dissipation at the boundary layer of the perfectly 
rigid crust and fluid core is given by [l| 



1 2 m+3 / 2 (m + l)! 2ttR 2 cPc [ Rc p(r) ( r\ 2m+2 dr 
Tv ~ 2fim(2m + l)!!2 m y r] c J 77 \R~c) ~R~c 

where the quantities R c , p c and r\ c are the radius, density and viscosity of the fluid at the 
outer edge of the core. In deriving expression ([5]) it is assumed that the crust is rigid and 
hence static in the rotating frame. The motion of the crust due to the mechanical coupling 
to the core effectively increases t v by a factor of (Av/v )~ 2 , where Av/v denote the difference 
between the velocities in the inner edge of the crust and outer edge of the core divided by 
the velocity of the core j^]. 

In neutron stars colder than about 10 9 K the shear viscosity is expected to be dominated 
by electron-electron scattering. The viscosity associated with this process is given by [l| 

Vee = 6.0 x 10 6 p 2 T- 2 , (g cm" 1 a" 1 ), (6) 

where all quantities are given in cgs units and T is measured in K. For temperature above 
10 9 K, neutron-neutron scattering provides the dominant dissipation mechanism. In this 



range the viscosity is given by [l| 

Vnn = 347p 9 / 4 T" 2 , (gcm^s- 1 ). (7) 

In the present work we neglect the effects of bulk viscosity, which are unimportant for 
T < 10 10 K. 

The fiducial viscous time scale f v is defined as 

T . =f .(^(TY 



where Qq — V^Gj) and p = 3M/47rR 3 is the mean density of the star. 
The gravitational radiation time scale is given by jl[ Q] 

1 32nGn 2m+2 (m-l) 2m fm + 2\ 2m+2 [ R < ,,.,,„_, 



p(r)r 2m+2 dr : (9) 



t gr c 2m+3 [(2m + 1)!!] 2 \m 

while the fiducial gravitational radiation time scale tqr is defined as 

2m+2 

n 

The critical angular velocity Q c , above which the r-mode is unstable, is defined by the 
condition r GR = r v and is given, for m = 2, by 



TGR =T G r[^) ■ (10) 



i t. 
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For a given temperature T and mode m the equation for the critical angular velocity, that is 
l/r(O c ) = 0, is a polynomial of order m + 1 in f2 2 , and thus each mode has its own critical 
angular velocity |25] . However, only the smallest of these (the m = 2 r-mode here) represents 
the critical angular velocity of the star and we concentrate our study on this r-mode. 

Moreover, the maximum angular velocity Qk (Kepler angular velocity) for any star occurs 



when the material at the surface effectively orbits the star 25j. This velocity is nearly 
Qk — l^o- Thus, there is a critical temperature below which the gravitational radiation 



instability is completely suppressed by viscosity and is given by jl| 

Employing Eqs. ( 11 ip and (|12p the critical angular velocity is expressed in terms of T c , that 
is 

Q c _Q K (T c \ 2/n _2 (T c \ 2/n 
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Once the equation of state for the neutron star core and crust is fixed, then all the ingredients 
of the r-mode instability, that is the transition density p c , the radial dependence of the mass 
density p(r), and the global properties of the neutron star (mass, radius and core radius) are 
determined in a self-consistent way. 



III. THE NUCLEAR MODEL 



he model used here, which has already been presented and analyzed in previous papers 



27H32|. is designed to reproduce the results of the microscopic calculations of both nuclear 



and neutron-rich matter at zero temperature and can be extended to finite temperature 



27|, 



2914321] . The energy per baryon at T = 0, is given by 



E b (n,I) = ±EP F tf" [(1 + J)5 /3 + (1 _ ^5/3] + 1_ A 



+ 



B [§ - (I + s s )J*] u° 



! + §£' [f-(I + x 3 )/ 2 ] u- 1 



(14) 



3 
2 



2=1,2 



3 x ^ 



i=l,2 



a 



Ci — 8Zi 
5 

d — 8Zi 



AA 3 /((1 + /) M ) 1/3 



F 



tan 



1/3' 



'((l-I)uf 3 + , ((1 - I)uf 3 



tan 



F 
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F 



In Eq. ( El) , / is the asymmetry parameter (/ = (n n — n p )/n) and u = u/uq, with no denoting 
the equilibrium symmetric nuclear matter density, n Q = 0.16 fm~ 3 . The parameters A, B, a, 
C\, C2 and B' which appear in the description of symmetric nuclear matter are determined 
in order that E^in = n Q ,I = 0) = —16 MeV, no = 0.16 fm -3 , and the incompressibility 
is K = 240 MeV and have the values A = -46.65, B = 39.45, a = 1.663, C x = -83.84, 
C2 = 23 and B' = 0.3. The finite range parameters are Ai = 1.5kp and A2 = 3kp and k F 
is the Fermi momentum at the saturation point no- The baryon energy is written also as a 
function of the baryon density n and the proton fraction x (x = n p /n), that is Ej,(n,x), by 
replacing 7=1 — 2x. 

The additional parameters x , x 3 , Z 1; and Z 2 employed to determine the properties of 
asymmetric nuclear matter are treated as parameters constrained by empirical knowledge 



271 ] . The parameterizations used in the present model have only a modest microscopic 



foundation. Nonetheless, they have the merit of being able to closely approximate more 




FIG. 1: The energy per baryon of symmetric nuclear matter (a) and pure neutron matter (b), 
as a function of the baryon density n, of the present model (MDIM) in comparison with those 
originating from realistic calculations. More details for the models UV14+TNI, UV14+UVII and 
AV14+UVII in Ref. Q and for the models A18+UIX and A18+du+UIX* in Ref. 34]. 



physically motivated calculations as presented in Fig. 1. More precisely, in Fig. 1 we com- 
pare the energy per baryon (for symmetric nuclear matter (Fig. la) and pure neutron matter 
(Fig. lb)) calculated by the present schematic model refereed as momentum-dependent in- 
teraction model (MDIM), with those of existent, state of the art calculations by Wiringa et 



al. 



and Pandharipande et al. 34]. 



A. Symmetry energy 



The energy E b (n, I) can be expanded around / = as follows 



E b {n, I) = E b (n, I = 0) + E sym ^n)P + £ S2/m , 4 (n)/ 4 + • • • + E symt2k (n)I 2k + . . . , (15) 



where the coefficients of the expansion are given by the expression 

1 d 2k E b (n, I) 



2k(n) 



(16) 



7=0 



(2A;) ! dl 2k 

In (fl~5"]) . only even powers of / appear due to the fact that the strong interaction must be 
symmetric under exchange of neutrons with protons i.e. the contribution to the energy must 



be independent of the sign of the difference n n — n p . The nuclear symmetry energy E sym (n 
is defined as the coefficient of the quadratic term, that is 

1 d 2 E b (n,I) 



E svm (n) = E svm o(n) 



(17) 

1=0 



The slope of the symmetry energy L at nuclear saturation density no, which is correlated 
with the crust-core transition density n t in a neutron star, is defined as 

L = 3n dEg y m{n 



dn 



n=riQ 

By suitably choosing the parameters Xq, x%, Z±, and Z 2 , it is possible to obtain different 
forms for the density dependence of the symmetry energy E sym (n) as well as on the value of 
the slope parameter L. We take as a range of L 50 MeV < L < 110 MeV where the value of 
the symmetry energy at saturation density is fixed to be E sym {no) = 30 MeV. Actually, for 
each value of L the density dependence of the symmetry energy is adjusted so that the energy 
of pure neutron matter to be comparable with those of existing state-of-the-art calculations 



33. 



34)- 



Fig. 2(a) displays the behavior of the nuclear symmetry energy as a function of the 
ratio u = n/riQ for various values of the slope parameter L. The aim of the above simple 
parametrization is to reproduce the nuclear symmetry energy originating from more realistic 
microscopic calculations and also covers the possible range of the nuclear symmetry energy 
dependence on the density. 



B. Proton fraction 



The proton fraction function of the baryon density n) in /3-stable matter. In this 

case we have the processes 

n — > p + e~ + v e p + e~ — > n + z/ e , (19) 

which take place simultaneously. We assume that neutrinos generated in these reactions 
have left the system. This implies that 

(i = ji n - ii p = /i e . (20) 
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FIG. 2: The nuclear symmetry energy (a) and the proton fraction (b) as a function of the ratio 
/no for various values of the slope parameter L. 



u = n 



The demand for equilibrium leads to equation 

d 



— ^E b {n,x) + E e {n,x)^=0, (21) 



or 

'dE h \ fdE f 



dx ) \ dx 1 

Finally, considering that the chemical potential of the electron is given by the relation (rel- 
ativistic electrons) 

H e = J k 2 F c 2 + m? e c A ~ k Fe c = hc(3rr 2 nx) 1/3 , (23) 



then Eq. ( 122]) is written 



dEi 



, hc^nx) 1 '*. (24) 

ox 

Eq. ( 1241) determines the proton fraction of /3-stable matter. In Fig. 2(b) we plot the proton 
fraction calculated from equation ( 1241) as a function of the ratio u = n/nQ for various values 
of the slope parameter L. According to Fig. 2(b) for low values of L, the neutron star 
matter consists mainly from neutrons and only a very small fraction of protons. However 
the increase of L (and consequently the increase of the stiffness of EOS) leads to an increase 
of the proton fraction. The values of the proton fraction define the kind of the cooling 
process of a hot neutron star. More precisely, it is well known that the direct Urea process 
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can occur in neutron stars if the proton concentration exceeds the critical value x cr u = 0.11 
for neutron star matter with electrons and x cr u = 0.148 when electrons and muons coexist 
in neutron star matter 35| . 



C. Nuclear equation of state for /3-stable matter 

The total pressure P(n,x), in the core of a neutron star, is decomposed into baryon and 
electron contributions 

P(n,x) = P b (n,x) + P e (n,x), (25) 

where 

2 dE b (n,x) 

P b {n,x) = n — . (26) 

The electrons are considered as a non-interacting relativistic Fermi gas and their contribution 
to the total energy density e e (n,x) and pressure P e (n,x) reads 



e e (n,x) = (37r 2 xra) 4/3 , (27) 



P e (n,x) = — (3n 2 xn) f \ (28) 

Now the total energy density e to t and pressure P to t of charge neutral and chemically equilib- 
rium nuclear matter is 

^tot = e& + e e , (29) 

P tot = P b + P e . (30) 

From Eqs. ( J29l) and ( J301) we construct the equation of state in the form e = e(P). When 
the electrons energy is large enough (i.e. greater than the muon mass), it is energetically 
favorable for the electrons to convert to muons 

e~ — > n~ + Up + v e . (31) 

However, in the present work we will not include the muon case to the total equation of state 
since the muon contribution does not alter significantly the gross properties of the neutron 
stars. 
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D. The thermodynamical method 



The core-crust interface corresponds to the phase transition between nuclei and uniform 
nuclear matter. The uniform matter is nearly pure neutron matter, with a proton fraction 
of just a few percent determined by the condition of beta equilibrium. Weak interactions 



conserve both baryon number and charge 36], and from the first law of thermodynamics, at 
temperature T = we have 

du = -Pdv - jldq, (32) 

where u is the internal energy per baryon, P is the total pressure, v is the volume per baryon 
(v = 1/n where n is the baryon density) and q is the charge fraction (q = x — Y e where x and 
Y e are the proton and electron fractions in baryonic matter respectively). In /3-equilibrium 
the chemical potential fi is given by fi = /i n — /i p = /i e where \i n and /i e are the chemical 
potentials of the protons, neutrons and electrons respectively. The stability of the uniform 
phase requires that u(v,q) is a convex function (37)]. This condition leads to the following 
two constraints for the pressure and the chemical potential 

_ « _ (!*) (*) > o, (33) 



V ® V ) a V ®Q J 



IT) >°- ( 34 ) 



dq 

It is assumed that the total internal energy per baryon u(v, q) can be decomposed into baryon 
(En) and electron (E e ) contributions 

u(v,q) = E b (v,q) + E e (v,q). (35) 



The relative theory has been extensively presented in our recent publication 



We consider 



the condition of charge neutrality q = which requires that x = Y e . This is the case we will 



consider also in the present study. Hence, according to Ref. [38( the constraints ( |33|) and 
031]), after some algebra lead to the following constraint 

C( n ) = 2r dEb{n ' x) + n 2 gM^) _ ( d 2 E b (n lX ) \ 2 / d 2 E b (n,x) \ ~' > 
dn dri 2 \ dndx J \ dx 2 J 

For a given equation of state, the quantity C(n) is plotted as a function of the baryon density 
n and the equation C(n) = defines the transition density n t . 
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FIG. 3: The transition baryon density rit (a) and the transition pressure Pt (b) as a function of the 
slope parameter L. For more details about the linear fit see text. 

IV. RESULTS AND DISCUSSION 

TABLE I: The slope parameter L (in MeV), the transition density nt (in fm~ 3 ), the transition 
pressure Pt (in MeV fm -3 ), the total radius R (in Km), the core radius R c (in Km), the core mass 
Mc in Mq and the central pressure P c (in MeV fm -3 ) correspond to a neutron star with mass 
M = 1.4M n . 



L 


rh 


Pt 


R 


R c 


M c 


Pc 


65 


0.0781 


0.317 


11.885 


10.956 


1.383 


86.4 


72.5 


0.0730 


0.319 


12.725 


11.631 


1.379 


63.2 


80 


0.0693 


0.295 


13.041 


11.898 


1.378 


57.5 


95 


0.0587 


0.155 


13.490 


12.385 


1.386 


49.2 


no 


0.0456 


0.0188 


13.646 


12.805 


1.397 


43.5 



We employ a phenomenological model for the energy per baryon of the asymmetric nuclear 
matter having the advantage of an analytical form. By suitably choosing the parametrization 
of the model we obtain various forms for the density dependence of the energy per baryon of 
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TABLE II: The same as in Table 1 for M = 1.8M, 



L 


n t 


Pt 


R 


R c 


M c 


Pc 


65 


0.0781 


0.317 


10.757 


10.287 


1.793 


331 


72.5 


0.0730 


0.319 


11.965 


11.328 


1.788 


176 


80 


0.0693 


0.295 


12.253 


11.584 


1.788 


159 


95 


0.0587 


0.155 


12.706 


12.052 


1.792 


133.5 


110 


0.0456 


0.0188 


12.951 


12.444 


1.798 


117 



neutron matter (see Fig. 1(b)), the nuclear symmetry energy (see Fig. 2(a)) and the proton 
fraction (see Fig. 2(b)), and in total the neutron star core equation of state. Additionally, 
the present model is applied for the determination of the transition density n t between the 
crust and the core. In order to complete the equation of state which describes the neutron 



star matter for densities lower than the transition density n t (equation of state of 



we employed the relative equation of state of Feynman, Metropolis and Teller 23] and also 



the crust) 



of Baym, Pethick and Sutherland 



24|. 



In Fig. 3(a), we plot the transition density n t , which is a fundamental quantity in the 
present study, as a function of the slope parameter L. A linear relation is found of the form 

n t = 0.1238 - 0.0007L. (37) 

According to relation fl37|) an increase of L (most stiff equation of state) leads to a decrease 
of n t and consequently to the extension of the radius core R c . In Fig. 3(b), we plot also the 
dependence of the transition pressure P t on L. The transition pressure does not affect directly 
the time scales corresponding to the r-mode instability. However it influences significantly the 
crustal fraction of the moment of inertia AI/I [39]. The ratio AI/I is particularly interesting 
as it can be inferred from observations of pulsar glitches, the occasional disruptions of the 
otherwise extremely regular pulsations from magnetized, rotating neutron stars [40|]. In 
Fig. 4 we display the mass-radius relation for neutron stars for the selected equation of 
states. All of them predict maximum mass neutron stars even higher than 1.8M . The case 
where L = 50 MeV has been excluded since it leads to very low value for the maximum 
mass of a neutron star. To further illustrate the mass-radius relation we present in Table I 
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FIG. 4: The mass-radius relations for the selected EOSs 

and II the transition density n t , transition pressure P t , the total radius of the star R, the 
core radius R c , the core mass M c as well as the central pressure P c for various values of the 
parameter L for neutron star mass 1.4M and 1.8M respectively . 

rRc 

The integral I(R C ) = / e(r)r 6 dr, which is a basic ingredient of the r-mode energy E 

Jo 

given by Eq. fl4| and also of the time scales (Eqs. (JSJ) and (jUJ) respectively) is plotted in 
Fig. 5, as a function of the slope parameter L for the interval 72.5 MeV < L < 110 MeV. 
e(r) is the energy density function defined as e(r) = p(r)/c 2 . It is obvious that the values of 
I(R C ) are correlated almost linearly with L in the considered interval. The least-squares fit 
expressions are respectively 

I(R C ) = ( 0.7 + 0.084L) 10 8 , (M = 1.4M ) (38) 

I(R C ) = ( 1.24 + 0.0845L) 10 8 , (M = 1.8M ). (39) 

The fiducial gravitational radiation time scale tgr combining Eqs. (Q and f fTOj) and after 
some algebra takes the form 



(40) 
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FIG. 5: The values of the fundamental integral I(R C ) = / e(r)r dr as a function of the slope 

J o 

parameter L for the interval 72.5 MeV < L < 110 MeV. For more details about the linear fit see 
text. 

where R,r are given in Km, M in M Q and the energy density e(r) in MeV/fm 3 . We can 
proceed further by replacing Eq. fl38l) in (j40]). In this way the time scale f GR corresponds to 
a neutron star with mass 1.4 M Q given by the simple expression 

fG « = 2 - 70736 (l^) 9 0.7 + 0.084L - (41 > 
where now R is given in units of 10 Km and L in MeV. It is worth to point out that tgr, 
according to Eq. f l4Tl) depends directly on the slope parameter L, as well as indirectly on the 
values of the radius R. In any case, Eq. fj4Tj) exhibits the equation of state dependence of 
the time scale tqr in a quantitative way. 

As a first approximation we can consider the case where the mass density of the neutron 
star p{r) is uniform that is p(r) = p = 3M/AttR 3 . Actually, this case is unphysical because 



the energy density does not vanish on the surface and the speed of sound c s = a/ dP/ dp is 
infinite 4l). However, we consider this case since it has been extensively applied for r-mode 
instabilities studies. After some algebra f^fP rox i s written 

-appro, _ , nK ( R YY 10K nA7 lil V" 



T GR 



195 (lOK^J {—) {-W) ■ < 42 > 
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In Fig. 6, as well in Table III, we compare the values of tqr and Tq^ ox . Actually both are 
increasing functions of the slope parameter L, for low values of L, but for higher values show 
a saturation trend. The approximated values for gravitational time scales are lower between 
8% - 16% (for M = 1.4M Q ) and lower by 23% - 28% (for M = 1.8M ). Consequently, the 
mean density approximation, concerning the gravitational time scales works better for low 
neutron star mass values. 

TABLE III: The fiducial time scales (in sec) for a neutron star with mass M = 1AM & 



L 




-approx 
T GR 


T Vee 


-approx 
T v ee 




-approx 

T V n n 


65 


2.417 


2.127 


30.855 


35.061 


69.735 


79.242 


72.5 


3.484 


3.176 


32.254 


35.380 


73.526 


80.656 


80 


3.986 


3.637 


33.833 


37.076 


77.637 


85.076 


95 


4.615 


4.123 


41.125 


46.035 


96.377 


107.882 


110 


4.468 


3.748 


56.977 


67.936 


137.840 


164.353 



TABLE IV: The fiducial time scales (in sec) for a neutron star with mass M = 1.8M 



L 


TGR 


-approx 
T GR 


T Vee 


-approx 

T Vee 


T V-nn 


-approx 

T Vnn 


65 


0.512 


0.366 


35.520 


49.747 


80.279 


112.433 


72.5 


0.870 


0.668 


36.686 


47.792 


83.629 


108.946 


80 


0.992 


0.760 


38.415 


50.145 


88.150 


115.068 


95 


1.182 


0.890 


46.613 


61.889 


109.238 


145.036 


110 


1.241 


0.895 


64.628 


89.664 


156.349 


216.920 



The fiducial viscus time f v which results from Eqs. ([5]), ([6]), ([7]) and ([8]) after some algebra 
is written for the case of viscosity due to electron-electron and neutron-neutron scattering 
respectively 
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FIG. 6: The fiducial time scales tqr, f ee , f nn as well as Tq^™, fe£ vrox , fnW r ° x as a function of 
the slope parameter L for a neutron star with mass M = 1.4M (a) and M = 1.8M@ (b). In case 
M = 1.4M© the gravitational times scales are multiplied by a factor 5 and in case M = 1.8Mq are 
multiplied by a factor 20. 



1.9 • W 



R 



3/4 



1M, 







1/4 



/lKm\ 6 1 1 f Rc 



\Rc J 



p c e c JO 



e(r)r 6 dr, 



(44) 



lKmy v M 

where p c and e c are the mass density (in gr cm -3 ) and the energy density (in MeV fin -3 ) 
respectively at the core edge. It is found that a linear relation holds between the mass 
transition density p c and the slope parameter L (and consequently a similar relation between 
e c and L) with the form 



Pc 



( 2.112 - 0.0121L) x 10 14 gr • cm" 



2.112 - 0.0121L x 56.1837 MeV • fin" 



(45) 



(46) 



Combining Eqs. (143)) . ( 144)) ( 145)) and (1461) we found for a neutron star mass with M = 1.4M Q 
the simple expressions 



7V, 



13.3053 



/ R \ 3/4 /10Km\ 6 1.24 + 0.084L , w 
(lOK^J (^r) (2.112- 0.01211)^ 



/ i? \ 3/4 /10Km\ D 1.24 + 0.084L , 
( TnrTZZ ) -5- ,0110 n.on,,, (M = 1.4M ) 



r, 



31.09 



"" V 1(IK| " / V /•'< / (2.112 - 0.0121 



(47) 
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FIG. 7: Temperature dependence of the critical angular velocity ratio £7 c /S7o for a neutron star 
with mass M = 1.4M© (a) and M = 1.8M (b) constructed for the selected EOSs. 

where now R and R c are given in units of 10 Km. By employing the mean density approxi- 
mation the corresponding expressions are 



-approx 
ee 



55.1657 



<T „ = 129 (^) 9/4 



9/4 



R c \ t M 
lOKm ) (uL 



3/4 



gr • cm 



3\ 3/2 



Rc 



M 



3/4 



gr • cm 



Pc,14 

3\ 13/8 



(49) 



(50) 



lOKm/ \IM & ; v Pc,i4 
where p c i4 is given in 10 14 gr • cm -3 . In Fig. 6, as well in Table III, we compare the fiducial 
viscous time scales. Both time scales are an increasing function of the slope parameter L, 
that is the more stiff is the equation of state the less viscosity effects on the r-mode. The 
approximated values for the viscosity timescales are between 9% — 16% (for M = 1.4M ) 
and between 23% — 28% (for M = 1.8M Q ). Consequently, the mean density approximation, 
concerning the viscosity time scales, works better for low values of neutron stars. 

In Fig. 7 are displayed the r-mode instability windows for neutron stars with mass 1.4M Q 
and 1.8M for the selected equations of state as a function of the temperature. For low 
values of temperature (T < 10 9 K) we plot the ratio 

f G ^ 2/i vio^y /n 



a 



n 



(51) 
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FIG. 8: The critical temperature as a function of the slope parameter L (a) and the transition 
density nt (b) for a neutron star with mass M = 1.4M and M = 1.8M . 

while for T > 10 9 K we plot the ratio 

The most striking feature is the location of the ratio Q c /Qo in a narrow interval (mainly in 
the case of neutron star with mass 1.4M ). Actually, the ratio fl c /^o increases around 7% 
(for T < 10 9 K) and around 9% (for T > 10 9 K) with the lower values corresponding to the 
case of L = 65 MeV and the higher to the case L = 80 MeV. It is concluded that the values 
of the ratio saturate for L close to the value 80 MeV. In the case of a neutron star with mass 
M = 1.8M the ratio fi c /^o increases around 10% both for T < 10 9 K) and for T > 10 9 K 
with the lower values corresponding to the case of L = 65 MeV and the higher to the case 
L = 80 MeV. 

The critical temperature T c , defined by Eq. (TTTI) . is plotted in Fig. 8, as a function of the 
slope parameter L and the transition density n t as well for neutron star with mass 1.4M 
and 1.8M . T c increases around 35% (for 1.4M ) and 45% (for 1.8M ) where the maximum 
corresponding to equation of state with L = 80 MeV. 

In Fig. 9 we compare the r-mode instability window for the five selected equation of states 
with those of the observed neutron stars in low-mass x-ray binaries (LMXBs) for M = 1.4M 
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FIG. 9: The critical frequency temperature dependence for a neutron star with mass M = 1.4M 
(a) and M = 1.8M & (b) constructed for the selected EOSs. In addition, the location of the observed 



short-recurrence-time LMXBs 
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43] 



and M = 1.8M . Firstly, we found that the instability window drops by ~ 20 — 40 Hz when 
the mass is raised from M = 1.4M to M — 1.8M . Following the study of Wen et al. [18j 
we examine four cases of LMXBs that is the 4U 1608-522 at 620 Hz, 4U 1636-536 at 581 Hz, 
MXB 1658-298 at 567 Hz and EXO 0748-676 at 552 Hz 
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431 ]. The masses of the mentioned 



stars are not measured accurately but the core temperature T is derived from their observed 
accretion luminosity. It is obvious from Fig. 9 that for a M = 1.8M all the considered 
LMXBs lie inside instability window. The case is similar for a M = 1.4M where the only 



exception the EXO 0748-676. According to discussion of Ref. [18] and finding in Refs 44j, |45 1 
the LMXBs should be out of the instability window. Consequently, one can presume that 
either the LMXBs masses are even lower than M = 1.4M or the softer equation of state 
is more preferred. However, additional theoretical and observation work must be dedicated 
before a definite conclusion. 

Another issue worth examining is the connection between the observed neutron stars in 
low-mass x-ray binaries and the nuclear physics input via the equation of state. It is a very 
interesting problem, if it is possible to constraint the nuclear physics input (for example 
the slope parameter L) from the related observation data. In general, this a very complex 
problem since the nuclear equation of state affects in different ways the possibility of an 
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r-mode instability and consequently additional work is needed as well to illustrate further 
this point. 

It is well known that for high temperature T > 10 10 K, which characterize mainly a new- 
born neutron star, the bulk viscosity is the dominant dissipation mechanism. However such 
a kind of dissipation is not considered in the present work. Actually, at high temperature, in 
a self-consistent treatment one has to consider also temperature effects on the nuclear equa- 



tion of state. It is known 



46| that temperature influences not only the density dependence 



of the equation of state but in addition both the transition density and the proton fraction. 
Actually the present MDIM model can be extended to include also the temperature effect 
on equation of state [28l-l32|. Work along these lines is in progress. 



V. SUMMARY 



In the present work we consider the effects, on r-mode instability, due to the presence of a 
solid crust in a neutron star. By employing a phenomenological nuclear model we calculated 
the equation of state of /3-stable matter which characterizes the neutron star core and is used 
for the location of the transition density at the inner edge between the liquid core and the 
solid crust. The stiffness of the equation of state parameterized via the slope parameter L, 
was varied on the interval 50 MeV < L < 110 MeV. The gravitational and the viscous times 
scales depend directly on the parameter L as well as indirectly on the transition density n t . 
As a consequence the critical angular velocity as well as the critical temperature depend on 
the equation of state. We found also that the instability window drops by ~ 20 — 40 Hz 
when the mass of a neutron star is raised from M = 1.4M to M = 1.8M Q . Finally we 
compared the r-mode instability window for the five selected equation of states with those 
of the observed neutron stars in low- mass x-ray binaries (LMXBs) for M = 1.4M Q and 
M = 1.8M . According to our finding one can presume that either the LMXBs masses are 
even lower than M = 1.4M Q or a softer equation of state is more preferred. 
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